The Markov Switching Multi-fractal models as a new class of REM-like models in 

1-dimensional space. 
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We map the Markov Switching Muhi-fractal model (MSM) onto the Random Energy Model 
(REM). The MSM is, like the REM, an exactly solvable model in 1-d space with non-trivial cor- 
relation functions. According to our results, four different statistical physics phases are possible in 
random walks with multi-fractal behavior. We also introduce the continuous branching version of 
the model, calculate the moments and prove multiscaling behavior. Different phases have different 
multi-scaling properties. 
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I. INTRODUCTION 

The Random Energy Model (REM) introduced by Der- 
rida pjj-tSil is one of the fundamental models of modern 
physics. Originally derived as a mean- field version of 
spin-glass models, it has subsequently been applied to de- 
scribe some features of 2-d Liouville models 0,0 j as well 
as the properties of quenched disorder in d-dimensional 
space with a logarithmic correlation function of energy 
disorder. 

The logarithmic correlation is easy to organize both for 
models on a hierarchic tree and in 2-d. The REM struc- 
ture in the d-dimensional logarithmic correlation case has 
been proved in Q using Bramson's results @ . In , 
a REM-like model in 1-d space has been solved directly, 
using the generalization of Selberg integrals [llj . A map- 
ping has been used in [l2| to map the REM onto strings. 

A REM can be formulated not only for the case of nor- 
mal distributions of energies, which corresponds to the 
logarithmic correlation function of energies on a hierar- 
chic tree, similar to logarithmic correlation functions in 
d-dimensional disorder case Q, but also for general dis- 
tribution of energies [4| . It is an open problem to find 
solvable models with the non-logarithmic correlation for 
energy disorder in 1-d space. To describe the fluctuations 
in financial market, [l4 | and (isj have constructed some 
dynamical models, the Markov switching multi-fractal 
models (MSM). The MSM has time translational sym- 
metry, contrary to cascade models, defined on hierarchic 
trees Q . The connection of the 1-d REM model [ll| with 
the multi- fractal random walk model [3- [13 was found 
in (igj . In this article we will prove that the dynamical 
models of [l3] and [l^ provide a 1-d REM where the 
correlation function for energies (In ut in our case) has a 
general character instead of being logarithmic. 

Let us give the definition of the MSM model, follow- 
ing [13, El, 113 ■ In the MSM model, one considers the 
sequence of variables rt, where t >0 describes a discrete 



moment of time: 

n = xtut, 
Xt has a normal distribution. 



(1) 



(2) 



and Ut is defined at the moment t of time via a product 
of k components M{t, I) 



Ut 



^Y[M{t,l). 



(3) 



1=1 



The variables M{t, I) are random variables with some 
distribution. 

Every moment of time our random variables are re- 
placed with new ones with a probability 



7, = 1 - cxp[-ab'''-% 



(4) 



where a > 0,b > 1,1 < I < k are parameters of the 
model. The parameter b plays the role of the branch- 
ing number in cascade models (models of random vari- 
ables on the branches of hierarchic models) and k is the 
maximal number of hierarchy on the tree. An impor- 
tant difference is that now 6 is a real number, while in 
case of hierarchic trees b should be an integer. Later we 
will formulate the continuous branching version of the 
model with a single relevant parameter V defined from 
the equation 



b^ 1 

L — > oo 



(5) 



We will use the notation L in sections II-H,III-B while 
investigating the multi-scaling properties of the model. 

The model is named as " random walk" model, because, 
due do Eq.(l), it is equivalent to the random walks with 
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an amplitude J^, when the time itself is a random vari- 
able see [13 and [l^ for a simple proof. The random 
variables M{t,l) are described via a Markov process, as 
any time period the transition probability depends only 
on the current state. There is a switching according 
to Eq.(4), thats why the authors of [l^ define it as a 
"switching" model. 

The distribution oi M{t,l) is chosen to ensure the con- 
straint 



< M{t,l) >= 1, 



(6) 



where <> means an average. 

We can take the log-normal distribution for M (i, /) or 
normal distribution for e^, defined as 



M(i,Z) = exp(/34), 



(7) 



where /? is similar to the inverse temperature in statistical 
physics. 

We consider a distribution for e: 



Vk He-Xf 
A = -/3V2. 



(8) 



We have for \t — t'\ < the following two-point correla- 
tion function for Ut : 



< UtUf >^ e 



(9) 



The correlation function for Inut, I nitf is logarithmic, as 
in models discussed in [loj . [§| . . [T6| . 

The previous expression is derived by observing that 
Ut and Ut' have identical M{t, I) for In \t — t'\/ In 6 levels 
of hierarchy. The probability that M{t, I) and M{f, I) 
are identical is 



exp[— |t — t'\e 



(10) 



Thus M{t,l) and M{t',l) are identical for 1-th level of 
hierarchy defined through the inequality 



k-'^<l<k. 
In 



(11) 



For the rest of the hierarchy levels M{t,l) and M{t',l) 
are different. 
When 



\n\t-t'\ 
V 



« 1, 



(12) 



the majority of the hierarchies have the same M(t, I) and 
M{t',l). 

Alternatively, when 



l-i^«l, 



(13) 



M{t, I) ^ M{t', I) for the the majority of 1. 

In section III we derive some more rigorous results. 



II. THE STATISTICAL PHYSICS OF MSM 
A. MSM with general distribution 

Let us consider a general distribution for e 

Po{q,£) = 7r- dhexp[-he + q(l){h)], (14) 

where g is a parameter indicating some effective length. 
We choose the distribution p(e) with some shift 



V 

= Po(-^,e- A) 



(15) 



to ensure the constraint given by Eq. (6): 

/oo 
depie) = e^^+^^(^) = 1. (16) 
-oo 



Thus we take 



A = - 



(3k 



(17) 



where k is our parameter describing the maximal hier- 
archy level. Hence, we have the following expression for 
the correlation function: 



< UtUf >^ e 



(y-ln|t-t'|/ln &)0(2/3)+20(/3)ln \t-t'\/\nb 



(18) 



B. The statistical physics versus the dynamics 



Let us define a partition function 

i—io l — l 



(19) 



where Al{i,l) is chosen from the distribution given by 
Eq.(15). 

Considering Eq.(l) as a dynamic process for a large 
period of time M, we define the probability distribution 



(20) 



We can define a statistical physics as well by considering 
z as a partition function for the 1-d model with quenched 
disorder. 

The average free energy, denoted as <lnZ > is: 



< InZ >=< lnz(io,e^) > 



(21) 



Let us consider a related model with standard distribu- 
tion for e given by po(e)j without the constraint of Eq.(4). 
We define 



i—iQ l — l 



(22) 
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where are defined through the distribution po(e) and 
the corresponding free energy is 

< InZo >=< lnzo(io,e^) > . (23) 

Eqs.(22) and (23) define a statistical physics model with 
configurations and special quenched disorder in 1-d 
space. Later we will focus on < InZ >. 
It is easy to check that 



< InZ >=< InZo > -V(j){P). 



(24) 



Eq.(24) is an exact relation, correct for any value of /3. 

It is easier to solve the model for < In Zq > . To cal- 
culate < InZo >, we will map the model onto the REM, 
and use the standard methods of REM. One can easily 
identify the most interesting transition in REM, from the 
high temperature phase to the SG phase, by just look- 
ing at the point in the high temperature phase where the 
entropy disappears. 

We will calculate the partition function's moments < 
Zq > and identify them with the < Z^^j^f >. 



where —Ei have independent distributions by given Eq. 
(14) with q = V, and Xi have normal distribution with 
variance 1. 

The moments of the partition function for this model 
can be calculated exactly by following [3] and 12j. 

These moments are identical to the expressions given 
by Eqs. (27) and (29). We assume that two models with 
identical integer moments have an identical free energy 
as well. 

The free energy of the REM model by Eq.(30) can be 
calculated rigorously following (l^ . 

At high temperatures, we have the Fisher zeros (FZ) 
phase with 

<lnZo>='-HZor = V^tm±l, (31) 

The transition point is at the point where the entropy 
disappears: 



0(2/3c) + 1 



(32) 



C. The moments in MSM model 

First of all we calculate 

The cross terms vanishes due to integration by xi. 
Let us consider now 



(25) 



<(Zo)2">=^..^<7.t,..ut,„ >. (26) 

While calculating the n-fold sum, we consider two prin- 
cipal contributions. 

The first case is when all t are close to each other and 
Eq.(12) is valid. There are A^i ~ such terms and the 
sum gives 

ln< (Zo)^" >=F + F0(2n^)+O(l). (27) 

The second case corresponds to the integration from the 
regions where U^i—ti are of the same order, and therefore 
the condition given by Eq.(13) is satisfied. There are N2 
such terms. 



In A^2 - nV 



V 



< 1. 



(28) 



As the vast majority of AP are different, the average gives 
In < (Zo)^" nV^ + nV^0(2/3) + 0(1). (29) 

D. The corresponding REM 

Consider now energy levels Ei and define the par- 
tition function 



Zrem — ^ • 



(30) 



Below this temperature the system is in the SG phase 
with the free energy 



W(2/3c). 



(33) 



Thus, we have found two phases. The phase given by 
Eq. (31) corresponds to the Fisher zero's phase, while 
that given by Eq.(33) is the SG phase. 



E. Asymmetric distribution 

So far we have considered the case of a symmetric dis- 
tribution of Xi . Let us now consider the asymmetric case 
described through the parameter 7, where 



< >= e 



(34) 



Now it is possible for the existence of a paramagnetic 
(PM) phase with the free energy 

< \nZo >= In < Zo >= (-7 + 1 + <P{P))V. (35) 



F. Large event 

Let us assume that at the starting moment of time 
there is a large event described through the parameter 
A, while for the other times Eq.(2) is valid. We consider 
the following partition function 

i—io + l l—l 

Now we can have the fourth, ferromagnetic (FM), phase 
with 



< InlZol >= In < Zq AV. 



(37) 
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Figure 1: The phase structure of the model with asymmet- 
ric distribution of weights. The case of normal distribution 



H. The scale dependence of the free energy 

Consider again the average distribution of Z, except 
that, instead of considering the sum over terms as in 
Eq.(20), consider the sum over I = e"^ terms, 



M/l 



(39) 



where e have the distribution by Eq.(14). At high tem- 
peratures we have Fisher zeros (FZ) phase with 



< hiZo >= V 
and for the critical point: 



/3c0'(2/?e) 



(40) 



(41) 



As a < 1, the /3c decreases with the decrease of a. 
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III. THE CASE OF CONTINUOUS 
BRANCHING 
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Figure 2: The phase structure of the model with a large 
event. The case of normal distribution (/f>(/3) — P^. 



Actually we can consider an infinite series of time, when 
after there is a member = — e"^^, while for other 
moments of time we calculate rt according to Eqs.(l) and 
(3). 



G. Transition points 

We should choose the proper phase by comparing the 
expressions given by Eqs. (31), (33), (35), (37) and then 
selecting the one which gives the maximum. 

For example, the system transforms from the FZ phase 
to the PM phase at 



,(l-7+0(/3))V ^ j2g-|: (1+0(2/3)) 



(38) 



where is the variance of Xt given by Eq.(2). 

In Figs. 1,2 we compare the numerics with our analyt- 
ical results for the free energy. 



A. The calculation of moments. 

All the formulae in the previous sections have been 
derived for the case of general values of b. Consider the 
case: 

V 



a = l,k = — , 6 = 1 + Sv, 
ov 



(42) 



and the random variables are distributed according 
Poi5v,e). 

The the i-th level of hierarchy is unchanged during the 
period of time t with a probability: 

expi-te"-^], 

V — idv. (43) 

Multiplying the probabilities in Eq.(43) for different lev- 
els of hierarchy, we obtain 

{[[{1 - exp(-te"))e^''2*(^) + exp(-te")e*''^(2;3)] ^ 
1=1 

exp[^ln((l - exp(-te''))e''''2'^('') + exp(-te'')e'''"^(2^)](]44) 

i 

Replacing the product by an integral and introducing 
variables Xi — ti/e^ , we derive 

dtidt2 < Ut^Ut^ >= 

^V{2{i+4>m) j dxidx2exp[[ dv[<i>2{e'',xi,X2)] 
Jo Jo 

$2(e^ xi, X2) = e-l^^-"^l^" (0(2/?) - 20(/3)) (45) 
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We get an asymptotic expression with the e ^ accu- 
racy in the Umit V ^ oo: 



re 

Jq dtidt2 < Utj^ut2 > 



^2V+2V<P{I3) 



(46) 



Similarly, we derive the expression for the multiple cor- 
relations: 



r,nV-\-n 



da;i..da;„e^o°°'i^*"(e":^i-^"). 



(47) 



The latter expression is 0(1), as has been assumed before 
in Eq. (29). 

For the 3-point correlation function we obtain 



p-(a:i2+2:23)j' 



(0(3/3) - 3m) 
+e-^i^^(l - e-^^«)(<?!)(2/3) - 2(/)(/3)) 
+ (1 - e-"i^^)e-"^^^(0(2/3) - 20(/3))] 



(48) 



where we denote Xi2 — \xi — X2\,X23 — \x2 — x^l- For n- 
point correlation function we need to consider 2"~^ terms 
in the expression of $„. 

We can identified this terms with different paths on 
a tree with branching number 2, the jumps to the right 
give a coefficient F{x, 1) and F{x, —1) for the left jump: 

F{x,l) = exp[-a;e"],F(a;,-l) = (1 -exp[-xe^]) (49) 

The path is fractured into clusters, when we have 1 sub- 
sequent right jumps. We define 



fi = m) 



(50) 



We should consider all the paths, the identify the n clus- 
ters of the given path with the length 1^ for the m-th 
cluster. Then we calculate 



$„(e",a;i...a;„) 



n-l 

^[nn^M+i,a)](E/'. 

paths i—0 m 



(51) 



B. Multi-scaling 



If we consider the model for / 



period of time and 



a normal distribution, we have in the high temperature 
phase 



F : 



< In Zo > _ 1 
aV ~ 2 



a 



(52) 




Figure 3: Free energy F versus an inverse temperature /3 
for the continuous branching model given by Eq.(42) and a 
normal distribution. The analytical results by Eqs.(52),(53), 
1^ = 14, A'^ = 100000 are given by smooth lines. 



and in the SG phase 

F^^^^^VVc 
aV 

The transition point is at 



(53) 



(54) 



In Fig. 3 we compare the numerics with our analytical 
results for the free energy. 

Let us calculate the multi-fractal properties of the 
MSM. We need to calculate the moments of the parti- 
tion function 



dti...dtn < ut^...ut^e 



where £,{n, (5) defines the multi-scaling, A{L) is some 
large numbers while C„ ^ 0(1). 

While calculating the moments, we slightly modify the 
formulae of the previous sections. Instead of Eq.(46) now 
we obtain 



/ 



dti..dtnU[ti)...u{tn) — I dxi-.dXn X 



I 



exp[ny(l + 0(/3)) / dw$„(-e",a;i...x„)]. (56) 

Jo ^ 

We consider the case \ L. Then, using the equa- 
tion 



$„(0, xi...Xn) = n0(/3) - (l){nP), 
we derive, integrating by parts: 



(57) 



f°° dy I 
/ — $„(y-,xi...x„) 
Ji y L 



= (</)(n/3) - n0(/3)) In ;^ - dy\ny<^'^{y,xi...Xn).m 
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Thus we get a multi-scaling with 

^{n,l3) =n + n4>{f3)-(t){n(3). 
Considering the moments of 



I 

t=i 



(59) 



(60) 



Consider again the model with z = exp[Q;T^] configura- 
tions. For the case 



we have, rescaling the result of [1), 

In < z^{a) >= nPV^/a 



(65) 



(66) 



where for ut we use the distribution given by Eq.(15), we 
obtain 



«(",/3)1"-e: / dxi...da;„e--''o°°'*'''"2^*"(^'^i-^'t^61) 



< 2;" > 

and here xi...Xn are time ordered. 

C. The moments for the model with random 
Boltzman weights. 

Let us calculate now the moments < z" > for 
I 



z = y^xtut- 

t=i 



(62) 



Using Eq.(8) from we derive 



(Ljy 



n\ I dxi...(ix„exp[- / In y$'„(y, 2/3, xi...x„)],(63) 



where there is a time ordering ti < t2 < --tn- 



In case of the multiscaling the right hand side is propor- 
tional to a. 

Thus there is a lack of any scaling in the SG phase and 
we can distinguish the SG and FZ phases checking the 
multiscaling property. We can distinguish FZ and SG 
phases also looking the tails of the distributions. 

For the FZ phase a simple rescaling of the results of 
111 gives for the large z: 



P{z) 



1 



Z ^ 



(67) 



while in the SG phase we used the rescaled result by Q 

1 



P{z) 



.1+ 



(68) 



where /3c/"\/2 is the transition point at a = 1. Eq.(68) is 
the result for the REM. For the logarithmic REM a more 
accurate expression includes a multiply In Z on the right 
hand side of Eq.(68) [22]. 



IV. THE DYNAMIC MODEL 



D. The multiscaling properties of different phases 

There are no simple order parameters to distinguish 
the FZ and SG phases. If we enlarge the free energy 
expression to the complex temperatures /3 = /3i -I- «/32, 
then in FZ phase there is a finite density p of partition 
function zeros, defined trough the formula [23| : 



p(/3i,/32) 



1 <lnz(/3i,/32)> 



27r 



V 



while in the SG phase this density is zero. 

SG and FZ phases have different schemes of replica 
symmetry (breaking) [2^ . 

There is a slow relaxation in THE SG phase. Unfor- 
tunately there are no any results about the dynamics of 
FZ phase to compare. 

More interesting is to distinguish two phases looking 
the multiscaling properties. We investigated well the 
muultiscaling properties of FZ phase. Let us investigate 
now the SG phase. 

Here there are the results by Gardner and Derrida Q 
about the the moments of partition functions and the 
probability distribution. 



Let us return to the dynamic model given by Eq.(l). 
Mapping 



-7^ 



(69) 



we identify this version of the model at /3 = as a finite 
time version of driven Brownian motion [25?]. In case of 
simple random walks with /3 = 0, there is a single phase. 
In the model given by Eq.(l) with /? ^ we have an 



(64) intrinsic large parameter L 



— 



describing the effec- 



tive number of configurations in the model. In this way 
the statistical physics enters into the dynamical prob- 
lem. Contrary to the driven Brownian motion and the 
Heston model [1^ , we have 4 different phases in the MSM 
model. The same is the situation with other models of 
multi-scaling random walks [T6| . 

To identify the choice amongst the two phases (FZ ver- 
sus PM) we consider 



(70) 



When C <C 1, the system is in the PM phase. Otherwise 
at C ^ 1 we have the FZ phase. 
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V. CONCLUSION 

We considered the dynamic Markov switching muhi- 
fractal models (MSM) and connected them with a new 
class of solvable statistical physics models of quenched 
disorder in one dimension. In these models there is both 
translational invariance and general distribution of dis- 
order. While the multiscaling in FZ phase has been in- 
vestigated before, nor the statistical physics properties 
properties, neither the phase structure has been inves- 
tigated before. We found the exact phase structure of 
the model. At different phases there should be different 
distributions P{z). In case of a symmetric random walk 
there are two phases in the considered model. At small 
parameters /3, the model is in the phase with non-zero 
density of Fisher's zeros. At high /3, the system is in 
the spin-glass phase, a pathologic phase with a slow re- 
laxation dynamics. For an asymmetric distribution of Xi 
there is a possibility for the third, paramagnetic, phase. 



Slight modification of the model allows the existence of 
the fourth, ferromagnetic, phase. It is possible to distin- 
guish different phases measuring the multi-scaling prop- 
erties of the model. The multi-scaling is broken in the 
SG phase. We also introduced a continuous hierarchy 
branching version of the MSM, gave expressions for the 
moments of the partition function, and calculated the 
multi-scaling indices. 

For applications it is important to calculate the frac- 
tional moments of the partition function. Perhaps we 
can use expressions for integer moments and use some 
approximate methods of extrapolation. Another inter- 
esting open problem is to investigate the dynamics of the 
model, looking for a new phase transition point in the 
dynamics, as is the case of the spherical spin-glass model 

DBS thanks Academia Sinica for financial support, 
Y. Fyodorov, S. Jain, Th. Lux, D. Sornette and O. 
Rozanova for discussions. 
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